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Abstract. - The symmetry of the superconducting order parameter in Bechgaard salts is still 
unknown, though the triplet pairing is well established by NMR data and large upper critical 
field H C 2(0) ~ 5 Tesla for H \\ a and H \\ b' . Here we examine the upper critical field of a 
few candidate superconductors within the standard formalism. The present analysis suggests 
strongly chiral f-wave superconductor somewhat similar to the one in S^RuCU is the most 
likely candidate. 



Introduction. - The Bechgaard salts (TMTSF) 2 PFg is the first organic superconductor 
discovered in 1980 [1]. For a long time the superconductivity is believed to be conventional 
s-wave [2], Recently the symmetry of the superconducting energy gap becomes the central 
issue [3,4]. The upper critical field at T = OK, _ff c2 ~ 5 Tesla for both H \\ a and H || b' 
for both (TMTSF) 2 PF 6 [5] and (TMTSF) 2 C10 4 [6] are clearly beyond the Pauli limit [7,8] 
indicating the triplet pairing. More recently the NMR data from (TMTSF) 2 PF 6 [9] indicates 
clearly the triplet pairing. Therefore the candidates for the superconductivity in Bechgaard 
salts are more likely within p-wave and f-wave superconductors. In the following we shall ex- 
amine the upper critical field of these superconductors following the standard method initiated 
by Gor'kov [10] and extended by Luk'yanchuk and Minccv [11] for unconventional supercon- 
ductors. Also we take the quasiparticle energy in the normal state as in the standard model 
for Bechgaard salts [2] 

£(fc) = w(|fc | — hp) — 2tb cos bk — 2t c cos ck (1) 

with v : vt, : v c ~ 1 : 1/10 : 1/300 and vt, — V2ti,b and v c = \^2t c c. There are earlier analysis 
of H C 2 of Bechgaard salts starting from the one dimensional models [12, 13]. However, those 
models predict diverging H C 2{T) for T — ► OK or the reentrance behaviour, which have not 
been observed in the experiments [5,6]. Also, the quasilinear T dependence of H C 2(T) in both 
(TMTSF) 2 PF 6 and (TMTSF) 2 C10 4 is very unusual. Among the models we have considered, 
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the chiral f-wave superconductor with A [kj ~ y^sgn{k a ) + i sinx 2 J cos % 2 , looks most 

promising, where xi =bk and \2 = ck where b and c are crystal vectors. 

Also if the superconductor belongs to one of the nodal superconductors [14] and if nodes 
lay parallel to k c within the two sheets of the Fermi surface, the angle dependent nuclear spin 
relaxation rate T^ 1 in a magnetic field rotated within the b' — c* plane will tell the nodal 
directions. 

Before proceeding, we show |A(fc)| os two chiral f-wave superconductors in Fig.^ a ) an d b). 

where |A(fc)| ~ [(1 + cos2 X i)(l - \ cos2 X2 )] * and |A(fc)| ~ [(1 + cos2x 2 )(l - \ cos2 X2 )] * 
for chiral fi and chiral f 2 respectively. 

Upper critical field for H \\ b' . - In the following we neglect the spin component of A(fc). 
Most likely the equal spin pairing is realised in Bechgaard salts as in Sr 2 Ru04 [4]. In this 
case the spin component is characterised by a unit vector d. Also d is most likely oriented 
parallel to c*. Let's assume d \\ c*, though H c2 {T) is independent of d as long as the spin 
orbit interaction is negligible. Experimental data from both UPt3 and Sr 2 RuC>4 indicate that 
the spin-orbit interactions in these systems are not negligible but extremely small [15]. We 
consider a variety of triplet superconductors (see Fig.^J: 

A. Simple p-wave SC: A(fc) ~ sgn(k a ). Following [16] the upper critical field is deter- 
mined by 

-lnt= r^-(l- Kl ) (2) 

J smh u 

-Clnt= / ^-(C-K 2 ) (3) 

J Q smh u 

where 



K 1 = ( e -P" 2 H 2 (1 + 2C P W)> (4) 




TMTSF) 2 PF 6 3 



Fig. 2 - Normalised H C 2(t) and C(t) for H || b' are shown in a) and b) respectively. Here solid, dashed 
and dotted lines are chiral f'-wave, chiral p-wave and simple p-wave respectively. 



K 2 = < e -"" 2 l s l 2 Qp 2 wV 4 + C (l - 8/m 2 | S | 2 + 12pV| S | 4 - yp 3 u 6 | S | 6 + ^ u 8 \s\ 8 ^) (5) 

and t — Tp-, p — VaV 2[2^T)^ i s ~ ~^2 S 9 n (ka) + *sinx2, X2 = ck and (. . .) means average over 
\2- Here v a , v c are the Fermi velocities parallel to the a axis and the c axis respectively. 
Here we assumed that A(r) is given by [16] 

A(r)~(l + C(a+) 4 )) (6) 

where) = £ C n e -eBx 2 -nk(x+tz)~ { ^- i s the Abrikosov state [17] anda+ = -J=(-id z - d x + 2ieHz) 
is the raising operator. 

Then in the vicinity of t ^ 1 we find p = ^^(-lnt) = 0.237697(- lnt) and C = 

93C(5) 
647C(3)^- 

For t — > on the other hand we obtain 

,. ,2 v a v c eH c2 (0) nl _ eq 

pn = hmpt = t— = 0.1583 (7) 

2(2nT c ) 2 

and C = —0.031. From these we obtain 

MO) = = 0-6659 (8) 

at l*=i 

Both /?o(i) and C(i) are evaluated numerically and shown in Fig. [21 a) and b) respectively. 
Here p (t) = t 2 p{t) = v v c eH c2 {t)/2(2irT c ) 2 . 

B. Chiral p-wave SC: A(fc) = 1/V2sgn(k a ) + zsin(x2)- Here -^=sgn(k a ) + i sin(x2) is the 

analogue of e 1 ^ in the 3D systems in the quasi ID system. For a chiral state the Abrikosov 
function is written as [18] 

A(r,fc)~( S + C S V) 2 )) (9) 
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where s = -^sgn(k a ) + isin(x2)- Then we obtained eq. |3]with 

K x = ( e - p " 2|s|2 (|s| 2 -2C|s| 4 )) (10) 
K 2 = < e -P" 2 N 2 (-| s | 4 + C \s\ 2 (1 - Apu 2 \s\ 2 + 2pV| S | 4 )) ) (11) 

and the same expressions for t, p,. . . 

For t -> 1 we find C = 1 - y/TE = -0.2247 and p = 0.3838(- hit). 
On the other hand, for t — ► we obtain C = —0.3660 and po = 0.27343. 
From these we obtain h(Q) — 0.71324. We obtain p(t) and C(t) numerically. They are 
shown in Fig. |21 a) and b) respectively. 

C. Chiral f-wave SC: A(k) ~ ds cos xi- H c2 (t) is determined from eq. |31where now: 

K x = ((1 + cos2 Xl ) e-"^ (\s\ 2 - 2pu 2 | S | 4 )) (12) 
K 2 = ((l + cos2 Xl )e-"" 2 l s l 2 (-pu 2 \ S \ 4 + C\ S \ 2 (l-4pu 2 \ S \ 2 + 2p 2 u 4 \ S \ 4 ))) (13) 

Here now (. . .) means the average over both xi an d Xi- Then it is easy to see that the 
chiral f-wave SC has the same H c2 {t) and C(t) as the chiral p-wave SC, since the variable Xi 
is readily integrated out. 

D. Chiral f-wave SC: A(fc) ~ ds cos xi- Now we have a set of equations similar to the 
chiral f-wave except 1 + cos2 X i in both eas. 1131 has to be replaced by | (1 + cos2x!). Then 
we obtain for t — > 1 C — -0.2247 and p = 0.5181(- hit). On the other hand, for t — > we 
find C = -0.3660 and p Q = 0.3734. 

We show po and C(t) of the chiral f-wave in FigJ2Ja) and b) respectively. 

Note that C(t) is the same for three chiral states (chiral p-wave, chiral f-wave and chiral 
f-wave) as well as chiral p-wave studied in [18] 

Therefore for the magnetic field H \\ b' , the chiral f-wave have the largest H c2 (t) if we 
assume T c and v, v c are the same. Also H c2 (t) of these states are closest to the observation. 

Upper critical field for H \\ a. - 

A. Simple p-wave SC: A (jzj — sgn{k a ). The equation for H c2 (t) is given by [16] and 
can be written as in eq|3|with 



K x = (e-^M 2 (1 + 2Cp 2 u 4 \s\ 4 )) (14) 
K 2 = ( e -P" 2 N 2 ( P V| S | 4 + c(l- 8pu 2 \s\ 2 + 12p 2 u 4 \s\ 4 - yA 6 M 6 + ^A 8 H 8 ))> (15) 

where t = p = '""^^ir)^ anc ^ s ~ sm Xi + «sinx2 with xi — bk and xi = cfc. 

Then for t -> 1, we find C = -^|yP and p = lni) = 0.2377(- hit). While for 

t -> C = Wo - \j{A) 2 + T2 = -0.0170129 and p = "° 2 "fff)| 0) = ^ exp [a + 2C(3 a ] = 

0.1751209, where a = —{In |s| 2 ) = 0.220051 and /3 = = £ - 1 = 0.0170. From these 
we obtain /i(0) = 0.73673. 

Both h(t) and C(t) are evaluated numerically and we show them in Fig. [31 a) and b) 
respectively. 
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Fig. 3 - Normalised H C 2(t) and C(t) for H \\ a are shown in a) and b) respectively. Here solid, 
dashed, dashed-dotted and dotted lines are chiral f'-wave, chiral f-wave, chiral p-wave and simple 
p-wave respectively. 



B. Chiral p-wave SC: A(fc) ~ y-^sgn{k a ) + zsin^j ■ Now H c2 (t) is determined by 
a similar set of equations as in sec. l.B. In particular we find for t — > 1 C = —0.027735 
and p = 0.212598(lni) while for t -> C = -0.067684 and p = 0.139672. We obtain 
h(0) = 0.6566. We show h(t) and C(t) in Fig. a) and b) respectively. 

C. Chiral f-wave SC: A(fc) ~ ds cos^i- Again we use a similar set of equations as those 
discussed in sec. l.C, we find for t ->• 1 C = -0.0356236 and p = 0.2 744495 (In t) while for 
t ^ C = 0.066 and p = 0.1920 and h(0) = 0.6997. Both h(t) and C(t) are evaluated 
numerically and shown in Fig.|3|a) and b). 

D. Chiral f'-wave SC: A(fc) ~ ds cos \2- Now we find for t — * 1 C = —0.05 and p = 
-0.2910(ln<), while for t -> C = -0.1019 and p = 0.2090. 

We have shown again h(t) and C(t) in Fig.|3|a) and b) respectively. 

Comparing these results with H c2 (T) from (TMTSF) 2 PF 6 and (TMTSF) 2 C10 4 [4,5], we 
can conclude both H || V and H || a the chiral f -wave SC is most consistent with experimental 
data. In particular these states have relatively large h(0) (see Table QJ. On the other hand 
almost the same H C 2(0) for H || b' and H \\ a has to be still accounted. 



Nodal lines in A(k). - We have seen that from the temperature dependence of H c2 (T), 
we deduce the chiral f-wave and chiral f are the most favourable. They have nodal lines on 
the Fermi surface (i.e. the Xi ~ Xi plane), the chiral f-wave SC at xi = if: while chiral 
f'-wave SC at X2 = ±f- 

These nodal lines may be detected if the nuclear spin relaxation rate is measured in a 
magnetic field rotated within the b' — c* plane. 

Following the standard procedure given in [14] the quasiparticle density of states in the 
vortex state for T << T c and E = is given by 
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Table I - Summary of results. Here po(0) = " 2 (2tt ffi an ^ ^(0) = ag^m * — 

— Si — 1* =1 

symmetry C(0) C(l) — gfjt = 1 po(0) ft(0) 

p-wave -0.031 0.2377 0.1583 0.6659 

H || 6' chiral p-wave -0.2247 -0.3660 0.3838 0.2734 0.71324 

chiral f'-wave -0.2247 -0.3660 0.5181 0.3734 0.72073 

p-wave -0.017 0.2377 0.1751 0,7366 

H || a chiral p-wave -0.066 -0.028 0.2126 0.1396 0,6566 

chiral f-wave -0.066 -0.035 0.2744 0.1920 0.6997 

chiral f-wave -0.1019 -0.05 0.2910 0.2090 0,7182 



where xw is the position of the nodal line on the xw axis. So for the chiral f-wave SC we 
find xio = 77 and N ^0, h) exhibits the simple angular dependence. On the other hand when 
nodal lines are on the \2 axis, the 9 dependence will be too small to see. Finally this gives 

Tfl (&) /T-J = (J^\ v 2 (eH) (1 + cos0 2 ) (17) 

for the chiral f-wave SC. 

We show the 9 dependence of T^ 1 in Fig. B for a few candidates. The chiral f-wave SC has 
the strongest 9 dependence (solid line) while the chiral h-wave SC (dashed line) and the chiral 
p-wave SC (dotted line) have a similar 9 dependence. 

Concluding remarks. We have computes the upper critical field of Bechgaard salts 
for a variety of model superconductors with the standard microscopic theory. We find: 
a) Assuming all these superconductors have the same T c , the chiral f'-wave SC (A(fc) ~ 
(-^sgn(K a ) + zsinx2^ cos X2) appears to be the most favourable with largest H^s for both 

H || b' and H \\ a; b) however, non of these states exhibit the quasi T linear dependence of 
H C 2(T) as observed in [4]; c) Also the present theory predicts H c2 (0) ~ (ww c ) _1 and (w6W c ) _1 
for H || b' and H \\ a respectively. This means H C 2(0) for H \\ a is about 5 time larger than 
the one for H \\ b' contrary to observation; d) from H C 2(0) ~ 5T and T c = 1.5K we can extract 
V 2 = s/v v c ~ 1.5 * 10 4 cms _1 , consistent with the known values of v, v c . 



* * * 



We thank Stuart Brown and Paul Chaikin for useful discussion on possible detection of the 
nodal structure of A(k)s in Bechgaard salts through NMR. We have benefited from discussion 
with Stephan Haas and David Parker. ADF acknowledges gratefully the financial support of 
Ministerio de Education y Ciencia (Spain) (AP2003-1383) and Jaime Ferrer and Paco Guinea 
for useful discussion. 



REFERENCES 

[1] Jerome D., Mazard A., Ribault M. and Bechgaard K., J. Phys (France) Lett, 41 (1980) 
L95 

[2] ISHIGURO T., Yamaji K. and Saito G., Organic Superconductors, Spring er-Verlag (Berlin 1999) 

[3] SlGRlST M. and Ueda K., Rev. Mod. Phys., 63 (1991) 239 

[4] Maki K., Haas S., Parker and Won H., Chinese J. Phys, 43 (2005) 532 



Ana Domi'nguez " JT '-~ TT ™™ — „™„ " yard salts (TMTSF)2PFg7 




e 

Fig. 4 - The angle dependent nuclear spin relaxation rate for a few nodal superconductors is shown. 
(Chiral f-wave, chiral h-wave and chiral p-wave are represented in solid, dashed and dotted lines.) 



[5] Lee I.J., Chaikin P.M. and Naughton M.J., Phys. Rev. B, 63 (2002) R180502 

[6] Oh J.I. and Naughton M.J., Phys. Rev. Lett, 92 (2004) 067001 

[7] Clogston A.M., Phys. Rev. Lett, 9 (1967) 266 

[8] Chandrasekhar B.S., Appl. Phys. Lett, 1 (1962) 7 

[9] Lee I. J., Brown S.E., Clark W.G., Strouse M.J., Naughton M.J., Kang W. and Chaikin 

P.M., Phys. Rev. Lett, 88 (2002) 017004 
[10] Gor'kov L.P., Soviet Phys. JETP, 10 (1960) 59 

[11] Luk'yanchuk I. and Mineev V.P., Soviet Phys. JETP, 66 (1987) 1168 
[12] Lebed A.G., JETP Lett, 44 (1986) 114 

[13] Depuis N., Mantambaux G. and Sa de Melo C.A.R., Phys. Rev. Lett, 70 (1993) 2613 
[14] Won H., Haas S. Parker D., Telang S., Vanyolos A. and Maki K., Lectures on the 

Physics of Highly Correlated Electron Systems IX, AIP Conference Proceedings 789 (Melville 

2005) 

[15] Maki K., Haas S., Parker D. and Won H., cond-mat/0504635, (2005) 

[16] Won H. and Maki K., Europhys. Lett, 30 (1995) 421 Phys. Rev. B, 53 (1996) 5927 

[17] Abrikosov A.A., Soviet Phys. JETP, 5 (1957) 1174 

[18] Wang G.F. and Maki K., Europhys. Lett, 45 (1999) 71 



